Polymer induced depletion potentials in polymer-colloid mixtures 
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The depletion interactions between two colloidal plates or between two colloidal spheres, induced 
by interacting polymers in a good solvent, are calculated theoretically and by computer simula- 
tions. A simple analytical theory is shown to be quantitatively accurate for case of two plates. A 
related depletion potential is derived for two spheres; it also agrees very well with direct computer 
simulations. Theories based on ideal polymers show important deviations with increasing polymer 
concentration: They overestimate the range of the depletion potential between two plates or two 
spheres at all densities, with the largest relative change occurring in the dilute regime. They under- 
estimate the well depth at contact for the case of two plates, but overestimate it for two spheres. 
Depletion potentials are also calculated using a coarse graining approach which represents the poly- 
mers as "soft colloids": good agreement is found in the dilute regime. Finally, the effect of the 
polymers on colloid-colloid osmotic virial coefficients is related to phase behavior of polymer-colloid 
mixtures. 
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I. INTRODUCTION 

Effective potentials are a key to unlocking the, equi- 
librium behavior of many soft-matter systemsSu. The 
basic philosophy behind this coarse-graining approach 
is that the initial effort in deriving these potentials is 
recouped when they are input into the well oiled ma- 
chinery of liquid state theoryQ, or when they are used 
in computer simulation^. An archetypal example is the 
depletion potential, induced between colloidal particles 
by non-adsorbing polymers. Asakura and OosawaO first 
showed that a bath of such polymers, characterized by 
their radius of gyration Rg, induces an attractive deple- 
tion interaction of range D w 2i?g between two plates. 
Their calculation was exact for non- interacting polymers. 
Later, the same authors, and independently VriJ^, de- 
rived a depletion potential between two colloidal hard 
spheres (HS) by approximating the (ideal) polymers as 
penetrable spheres. This is often termed the Asakura- 
Oosawa (AO) model. 

A good example of the effective potential coarse- 
graining approach is the calculation of the phase be- 
havior |-|0f polymer-colloid mixtures by|-Gast, Hall, and 
RusselL^nd also by Lekkerkerker et aVB and Meijer and 
FrenkelE2l. They found, using an AO depletion potential 
approach, that the fluid-fluid phase line of colloids of ra- 
dius Rc becomes metastable w.r.t. the fluid-solid phase 
line if the size ratio q — Rg/ Rc is less thajji_about 0.35, in 
qualitative agreement with experimentsH'EEI. Their work 
demonstrates how an accurate knowledge of the deple- 
tion potential can lead to a good understanding of the 
equilibrium behavior of colloid polymer mixtures. The 
latter are important not only because of their relevance 
to many industrial and biological processes, but also be- 
cause they form an important model system for equilib- 



rium and non-equilibrium behavior in soft matter science. 

Whereas the depletion interaction for ideal polymers 
is now quantitatively understood, the depletion interac- 
tion induced by polymers with excluded volume inter- 
actions is only qualitatiyely understood. . Experiments 
on the phase behavioiO'EEl and structuretlHlj of poly- 
mer colloid mixLutfis also show deviations from the sim- 
ple AO modellj^O, [-as-do direct measurements of the 
depletion potentialsEjilJ. Theoretical attempts to di- 
rectly calculate the depletion potentials for interacting 
polymeis include scaling theprylla, self-consistent fieW 
theoryO, pejvUirbation thcoryEi'EJ, direct simulationstj, 
RG theoryc3c3, as well an interesting new "overlap ap- 
proximation" methodO. All these approaches show sig- 
niflcant deviations from ideal polymer behavior, but 
many questions still remain. This is in contrast to bi- 
nary HS colloid mixtures, where the deviations from 
the AO potential can now be quantitatiiiely calculated 
with density functional theory (DFT)cJ Ea, and the ef- 
fects of nonjdeality on the phase behavior are fairly well 
understood^. The goal of our paper is to derive a theory 
of similar accuracy for the depletion potential induced by 
interacting polymers in a good solvent 

Before we proceed, an important caveat is that the de- 
pletion potential becomes less relevant for phase behav- 
ior at larger size ratios q, since-maxK-body interactions 
become increasingly importantE2rEjO. When q » 1, 
i.e. when the colloids are much smaller than the poly- 
mers, other approaches, which treat the polymers on a 
monomer level, are more|¥elevant. ExampleSiiHclude inte- 
gral equation techniquegSil, scaling theoryBEj, or renor- 
malization group theory (RG)EJ. In this paper we con- 
centrate on the regime where Rg is of the order of Rc or 
smaller. 

We have recently performed a systematic study of the 
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insertion free-energy ofpa,single colloidal particle in a bath 
of interacting polymersLj (henceforth referred to as paper 
I), and found important deviations from ideal polymer 
behavior. In particular, we found that the range of the 
depletion layer decreases with increasing density p, and 
that this effect is most pronounced in the dilute regime, 
where p < p*; here p* — |7ri?g is the overlap density. In 
the semi-dilute regime, where p/p* > 1, we found that 
the effects of thp interactions could be well described by 
scaling theoryEj. Since these one-body effects were not 
well captured by (non- additive) HS like models, we do 
not expect straightforward extensions of the DFT tech- 
niques that work so well for binary hard HS mixtures, 
to also perform well for two-body depletion potentials in 
polymer-colloid mixtures. 

To calculate the depletion potentials, we use a similar 
approach to that used in paper I - a combination of scal- 
ing theories and computer simulations. Since we found 
in paper I that the range of the depletion layer decreases 
with increasing density, we expect a similar trend for the 
range of the related depletion potential. Throughout 
work we focus on the dilute and semi-dilute regimea23 
of the polymers, where the monomer density c is low 
enough for detailed monomer- monomer correlations to 
be unimportant; the melt regime, where c becomes ap- 
preciable, will not be treated here. Since our models are 
all athermal, we set /3 = l/(fcsr) = 1. 

The paper is organized as follows: The case of the de- 
pletion interaction between two plates is discussed in de- 
tail in section II. We show that it is closely related to the 
problem of determining the surface tension 7^, (p) , which 
was solved in paper I. Just as was found for 7-u,(p), the 
depletion potential simplifies in the semi-dilute regime. 
In section III, we discuss the depletion interaction be- 
tween two spheres. We compare the results of direct 
Monte Carlo (MC) simulations of the interaction between 
two HS colloids induced by a bath of self avoiding walk 
(SAW) polymers to a_cotential derived within the Der- 
jaguin approximationE]. The Derjaguin approximation 
works much better than one would naively expect be- 
cause of a cancellation of errors related to the defor- 
mation of polymers around a sphere. Using an exten- 
sion of the Derjaguin approximation, we derive a new 
semi-empirical depletion potential which appears to be 
nearly quantitative for q < 1, the regime where deple- 
tion potentials are most relevant for the phase behavior 
and structure of polymer-colloid mixtures. We derive the 
scaling behavior with density, and find important devi- 
ations from the AO model and other ideal-polymer the- 
ories. We also calculate the depletion potential between 
spheres within OjUx new "polymers as soft colloids" coarse- 
graining schemeLj Cj, where each polymer is represented 
by a single particle, interacting via a density-dependent 
effective potential. Here we again find good agreement 
with the direct MC results for the dilute regime, but for 
p/p* > 2, significant deviations are found. In section IV, 
we discuss the effect of polymer density on the virial co- 
efficients between two colloids, and relate this to phase 



behavior of polymer-colloid mixtures. 



II. DEPLETION POTENTIAL BETWEEN TWO 
WALLS 

A. Surface tension near a single wall and the 
depletion potential at contact 

Immersing a single hard wall (or plate) into a bath 
of non-adsorbing polymers in a good solvent reduces the 
number of configurations available for the polymers. This 
in turn results in an entropically induced depletion layer 
p(z) near the wall, discussed in more detail, for exam- 
ple, in paper I. Associated with the creation of this de- 
pletion layer is an intcrfacial free energy cost per unit 
area A, i.e. a wall- fluid surface tension 7„(p), which typ- 
ically depends on the bulk|-density p, or equivalently the 
bulk chemical potential pEj. Bringing two such walls, 
of area A, together from an infinite distance apart to 
a distance x where the two depletion layers begin to 
overlap, changes the total free energy of the polymer 
solution. This change in the free energy or grand po- 
tential per unit area is called the depletion potential 
W{x) EE {n{x) - n{x = 00)) /AB. When x = 0, i.e. when 
the two walls are brought into contact, the depletion po- 
tential reduces to the simple form 



W{0) = -27^ (p). 



(1) 



reflecting the fact that the two depletion layers are com- 
pletely destroyed. 

In paper I, we used an extension of the Gibbs adsorp- 
tion equation to express the surface tension near a single 
wall as 



(^|r(.')V. 



(2) 



The cLesiMition of this equation can be found, for exam- 
ple, irc30. Here H(p) is the osmotic pressure of the poly- 
mer solution, and r(p) is the reduced adsorption near a 
single wall, defined as: 



p dp 



P(f) 
P 



-l\dz 



(3) 



where fi'^^/A is the excess grand potential per unit area. 
In paper I, we used computer simulations of a SAW poly- 
mer solution in a good solvent near a wall to calculate 
p(z) and r(p) for several values of p/p*. As discussed in 
paper I, the SAW on a cubic lattice is a very good model 
for polymers in a good solvent. In the scaling limit, where 
the length L tends towards 00, its properties are univer- 
sal, and agree with experiments on polymers in the same 
good solvent regimcHEj. For example, the radius of gy- 
ration scales as Rg ^ L'^ , where i' is the Flory exponent, 
taken to be = 0.588 in this paper. By using an accurate 
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fitting form which takes into account the correct scaling 
behavior, we expressed T(p) for all densities in the dilute 
and semi-dilute regimes. When this was combined with 
an RG expression for the pressure, we were able to cal- 
culate the surface tension "fw{p) as a function of density 
through Eq. (H). Our ijesults agreed very well with some 
recent RG calculation^. 

Here we use jw{p) together with Eq. (|^) to directly 
calculate the depletion potential at contact, as shown in 
Fig. |l|. These results are also compared to dijy?qt com- 
puter simulations of L = 100 SAW polymertOo. Even 
though the L — 100 simulations have not quite reached 
the scaling limit, the well-depths at contact are remark- 
ably well reprnduced, a result also found by Tuinier and 
LekkerkerkeinJ with a similar theory. 
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FIG. 1. Depletion potential between two walls or plates 
at contact VF(0) = —2'yw{p) per unit area Rg. The full lines 
result from Eq. (|^) , while the long-dashed lines come from the 
simpler semi-dilute scaling expression (^. The diamonds de- 
note previously published L = 100 SAW simulation resultsE3. 
The dashed line represents the simple ideal polymer form 
Wid{0), given by Eq. (^, while the dot-dashed line shows the 
naive improvement obtained by substituting the true polymer 
osmotic pressure for the ideal polymer pressure. These two 
Asakura-Oosawa type approximations bracket the correct re- 
sult. The inset shows the strength of the simpler semi-dilute 
scaling expression relative to the the full potential at contact 
(dotted line). The two coincide for higher densities. In the 
low density limit, the semi-dilute scaling expression overesti- 
mates the true value at contact by a factor 1.5. 



B. A simple theory for the depletion potential 

The simplest approximation for the depletion poten- 
tial at finite separation x would be a linear form with a 
slope equal to the osmotic pressure. This follows because 
n(p)x is the volume term proportional to the work per 
unit area produced by the osmotic pressure: 



W{x) = W{0) + n{p)x; x<D^ (p) 
W{x)=0 ; x>D^{p) (4) 



where the range is given by 
W{0) 



272 (p) 



n(p) n(p) 



(5) 



We note that this approximation is similar to that 
adopted by Joanpv, Leibler, and de Gcnnes who, in their 
pioneering papertS, approximated the force between two 
plates as constant for x < 7r^(/o) and zero for x > n£^{p), 
where ^(p) is the correlation length, the relevant length- 
scale in the semi-dilute regimetJCJ. This also results in 
a linear depletion potential. 

In Fig. H we compare our simple linear potential to the 
direct SAW simulations. The overall agreement is strik- 
ing. The only (small) deviation occurs at larger distances 
X where the true potential rounds ofLand develops a very 
small maximum before going to zercc3. The range (p) 
of this simple potential decreases with p/p*. In fact, as 
shown in Fig. ^, the largest relative rate of decrease oc- 
curs in the dilute regime, so that at p/p* — 1, the range 
\s D — 1.25Rg, about 58% of the low density limit of 
D w 2.15. In general, the rapge of a potential is not al- 
ways unambiguously definedc^l. Of course for our linear 
one it is, and this simple definition would seem a very 
reasonable definition for the range of the full depletion 
potentials depicted in Fig. 0. 




O pjp* 


= 0.26 


A p,/p* 


= 0.62 


O p,/p* 


= 1.01 



0.5 



x/R„ 



1.5 



FIG. 2. Depletion potential between two walls or plates 
W{x) per unit area Rg. The symbols denote simulations for 
L = 100 SAW polymers, the straight lines are the simple 
theory of Eq. (M), which provides a near quantitative fit to 
the simulation data. The range Dm(p = 0) « 2.15 is shown 
as a vertical line. Note how much the range decreases with 
density, even for these results in the dilute regime. 
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C. Scaling theory for the depletion potential in the 
semi-dilute regime 

Further simplifications occur in the semi-dilute regime. 
For example, when the scaling formsES for the osmotic 
pressure, 11 ^ p3i'/(3i'-i) ^t^^ reduced_adsorption 



are used in Eq. 



then, as 



shown in paper I, the integral simplifies and the poten- 
tial at contact takes on the very simple form: 

= 3n(p)f(p) ^ « (6) 

The linear depletion potential (^) then reduces to: 

Wsdix) = n{p) (sf (p) + ; x<Dsd 
Wsdix)^0 ; h>D,d, (7) 

and the range has simplified to 



(8) 



As can be seen in Figs |l| and ^, these simple expressions 
for the well-depth and the range work remarkably well in 
the semi-dilute regime, especially for p/p* > 2. 



D = —2T{p), so that Eq. also overestimates the range 
at low densities by a factor 1.5, as can be seen in Fig. y. 

In |the semi-dilute regime one can identify r(p) w 
— Therefore the ansatz D^Ju)) — 7rf(p), orig- 
inally postulated by Joanny et al.ll3, is very close to 
Dsd[p) = — 3r, the expression we derived for the range of 
the depletion potential in the semi-dilute regime. Their 
potential is therefore quite accurate in the semi-dilute 
regime, but overestimates the range and the well-depth 
by a factor slightly larger than 1.5 in the dilute limit. 



D. Comparison with theories for non-interacting 
polymers between two walls 

For completeness we compare the results obtained in 
the previous section to those for ideal polymers, first ob- 
tained by Asakura and Oosawa in 1954la. Their (exact) 
depletion potential can |be quite accurately approximated 
by a simple linear formtJ: 



W,d{x)=(} 



X > —=Rn 



(9) 



D = -2yjn 
D.. = -3r 
D for ideal polymers 




FIG. 3. Comparison of the range Dw{p) of the deple- 
tion potential with the simpler semi-dilute scaling expressions 
Dsd{p), both in units of Rg. As density increases well into the 
semi-dilute regime, the two expressions coincide. In the low 
density limit the semi-dilute expression is exactly 1.5 times 
the full expression. Note how much the range differs from 
the density independent result for non-interacting polymers 
Ad = ^ 2.26 



Deviations do occur for low densities since according 
to @) and (I), W{x = 0) « 2n(p)f (p) for p ^ 0, the 
same form as for ideal polymers. Eq. (ra) therefore over- 
estimates the well depth by a factor 1.5, as can be seen 
in the inset of Fig. |l]. In the same limit the range of 
the depletion potential reduces to the functional form 



However, for interacting polymers, this is only true in 
the limit p — > 0; the validity of this expression rapidly 
deteriorates with increasing density. As shown in Fig. 
Eq. (||) underestimates the well depth for all but the low- 
est densities, while, as shown in Fig ^, it overestimates 
the range at all densities. One might think that replacing 
the ideal pressure n(p) = p in Eq. @ by the pressure of 
an interacting polymer system would bring an improve- 
ment for the well depth, even if this does not improve 
the approximation for the range. Instead, as shown in 
Fig. y, this naive approach leads to a severe overestimate 
of the well-depth. In fact, for the semi-dilute regime, the 
contact value of this naive "improvement" would scale 

as WX''"=(0) = (4/v^)n(p) - p3''/(3'>-l) ~ p2.309 

stead of the correct W{0) ^ p2''/(3''-i) ~ pi-539 scaling. 
Even for p/p* = 1 the differences are significant: W{0) — 
-0.93, while Wid(O) = -0.54 and W^f'^^iO) = -1.68. 
For p/p* = 10, WiO) = -25.4, WidiO) = -5.4 and 
— —208! In other words, the simple heuris- 
tic ansatz is almost never a real improvement; for inter- 
acting polymers, the Asakura-Oosawa potential between 
two plates is only accurate at the very lowest of densities. 



III. DEPLETION INTERACTIONS BETWEEN 
SPHERES 

We now turn to the polymer induced depletion inter- 
action between two hard spheres. Similarly to the case 
of two walls, when the surfaces of two such spheres are 
brought to within a distance where their depletion lay- 
ers begin to overlap, the total free energy of the system 
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changes; this change again defines the depletion poten- 
tial. 5 
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FIG. 4. Depletion forces between two spheres for three size 
ratios q, as a function of x, the distance between the surface 
of the two spheres. The symbols denote L = 500 SAW com- 
puter simulations for three densities: p/p* = 0.58 (circles), 
p/p* = 1.16 (squares), p/p* = 2.32 (diamonds), the solid lines 
are to guide the eye. These are compared to results from the 
Derjaguin approximation for the same size ratios and densi- 
ties: p/p* = 0.58 (dotted lines), p/p* — 1.16 (dashed lines), 
and p/p* = 2.32 (dot-dashed lines). 



Since we found in paper I that the surface tensions as- 
sociated with the polymer depletion layers around a sin- 
gle sphere show behavior which differs from that of ideal 
polymers or that of the AO model, we expect to find qual- 
itative differences in the depletion potentials as well. The 
trends are expected to be similar to those foimd for the 
depletion potential between two walls, namely that the 
range should decrease and the well depth should increase 
with increasing polymer concentration. 



A. Full SAW polymer simulations 



To directly calculate the polymer induced depletion 
potential between two spheres, we performed grand- 
canonical simulations of L = 500 SAW polymers on a 
lattice of size 240x150x150, and computed the osmotic 
pressure or force exerted on two hard spheres placed in 
the same simulation box. The configuration space of 
polymers was sampled in the grand canonical ensemble 
with a combinaticin of Configurational Bias Monte Carlo 
and pivot movesu. The depletion force on the spheres 
was obtained from the ratio of the acceptance of virtual 
inward and outward moves of the spheres. 

In Fig. U these forces are shown for three different size 
ratios q — Rg/Rc, namely q = 0.67, q ~ 1.05, and 
q — 1.68. For the first two size ratios, we computed the 
forces at three densities, p/p* = 0.58, p/p* = 1.16, and 
p/p* = 2.32c3, for q — 1.68, this was only done for the 
lower two densities. As expected, for a fixed size ratio, 
the range decreases, and the force at contact increases 
with increasing polymer density. For a given density, the 
range appears to contract slightly with increasing q, as 
might be expected, since the polymers can deform more 
readily around the smaller colloids (see the Appendix for 
further discussion of this point). In each simulation we 
keep the size of the polymers fixed, so that the larger 
size ratios essentially correspond to smaller spheres. The 
computational costs scales with the size of the simulation 
box, which roughly sets the number of polymers needed 
to achieve a density p/p* in the accessible volume left by 
the colloids. 

The force can then be integrated to obtain the effective 
depletion potential: 



V{x) 



F{y)dy. 



(10) 



Results are shown in Fig. |g for the same parameters as 
in Fig. ^. Because the curves are integrated, they appear 
smoother than the forces. There may still be some resid- 
ual error in the potentials, which may explain why the 
range for the highest density seems slightly larger for the 
potential than for the force. Also, the errors are likely 
too large to decide whether or not there is a repulsive 
bump in tie potential. If it does exist, it is probably 
very smallcZl. This is in contrast to hard-core systems. 
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where such bumps can be pronouncedSilS. The weak- 
ness or absence of a repulsive barrier in the case of inter- 
acting polymers can probably be traced to the very low 
monomer concentration c. For shorter polymers, whose 
behavior deviates significantly from the L —^ oo scaling 
limit, and where a substantial monomer density c can 
more easily be achieved, such repulsive bumps might oc- 
cur more readilyEj. 
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FIG. 5. Depletion potentials between two spheres for the 
same set of parameters as in Fig. ^ These are compared to 
the semi-empirical depletion potential Vs{r) (Eq. 0) for the 
same size ratios and densities. 



As expected, the range decreases while the force at con- 
tact and the potential well-depth increase with increasing 
polymer density. For a fixed density and polymer size Rg , 
the force at contact and the well-depth decrease as the 
spheres become smaller (or q becomes larger). The same 
effect occurs for ideal polymersBEl, where it has a simple 
geometric origin: the volume of a depletion layer of a 



given width ^ 
hard spheres. 



Rg decreases with decreasing size of the 



B. The Derjaguin approximation for interacting 
polymers 

In the case of two walls, the well-depth at contact has a 
clear interpretation in terms of the complete destruction 
of two depletion layers, and can therefore be expressed 
as a simple function of the wall-polymer surface tension 
jw{p), as shown in Eq. (0). For two spheres, the deple- 
tion layers do not completely overlap, and the depletion 
potential at contact is related to the surface tension of 
polymers surrounding two spheres in a dumbbell type 
configuration. A simple expression for the well depth at 
contact in terms of the surface tension ^s{p) is therefore 
less obvious. 

One way to make contact with the two wall case is 
to use the Derjaguin approximation^, which relates the 
force between two spheres of radius Rc to the potential 
between two walls, W{x), in the following way: 



(11) 



where x is the distance between the surfaces of the two 
spheres. In principle this approximation should only be 
accurate for very small size ratios g, i.e. for very large col- 
loids. However, as shown in Appendix A, we expect there 
to be a cancellation of errors, related to the deformation 
of the polymers around spherical particles, which makes 
the Derjaguin approximation work much better than one 
would naively expect. This is confirmed in Fig. ^ where 
the Derjaguin expression for the force, taken from Eq. (^) 
and Eq. (pT|), is shown to be a surprisingly good approx- 
imation. It is most accurate for the smallest q, as ex- 
pected, but even for q = 1.68, where normally one would 
not expect the Derjaguin approximation to be useful at 
all, it is still reasonable. 

By combining Eq. (|) with Eqs. @ and (|ll|) we ob- 
tain the following Derjaguin expression for the depletion 
potential between two spheres 

y^-^- (x) ^ -^rMp) {d^ (p) -xf (12) 

for < x < D^{p)- V^''''^{x) = for a; > D^{p). By 
the nature of the Derjaguin approximation, Djj^{p) is the 
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same range as found for two walls in Eq. (^. This sim- 
ple potential is compared in Figs. || to the direct L = 500 
SAW simulation results for q = 1.05. The correspondence 
is surprisingly good given the large size ratio 5, although 
it is not quantitative as in the wall-wall case. Similar re- 
sults are found for the two other size ratios. Overall, the 
Derjaguin approximation overestimates the well-depth, a 
consequence of the slightly longer ranged forces found in 
Fig. 1^. Again, the cancellation of errors found for the sim- 
pler AO model in Appendix A helps explain why Eq. jl^ ) 
works reasonably well in a regime where the Derjaguin 
approximation would normally break down. 
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0.5 1 1.5 

FIG. 6. Comparison of the L = 500 SAW computer simu- 
lations for the depletion potential with the Derjaguin approx- 
imation of Eq. (hd) for q — 1.05. The symbols are the same 
as in Fig. M. 



C. The Derjaguin approximation in the semi-dilute 
regime 

We now turn to the scaling behavior of the depletion 
potential in the semi-dilute regime. By using the sim- 
ple expression for the range Dsd{p) given in Eq. (||), the 
Derjaguin approximation for the depletion potential (IH 
reduces to 



Derj 



sd 



{x) 



R,n{p) 3f(p)+x 



(13) 



for X < Dsd{p)- Vjd'^^'' i^) = for a; > Dsd{p)- For two 
walls, the simplified expression (H) for the range overes- 
timates the true range by a factor 1.5 in the low-density 
limit. Here the overestimate is slightly larger, since, as 
shown in the appendix, the deformation of the polymers 
around a single colloid reduces the depletion layer width 
for decreasing sphere size Re- Similarly the well-depth 
at contact is also overestimated in the low density limit, 
but now by a factor limp^o V;^'"'^ (0)/F^'=''^ (0) = 2.25. 
In the semi-dilute regime the two expressions come closer 
for increasing density. For example, at p/ p* = 2 the over- 
estimate at contact is only V;^'"'^ (0)/^^^''^ (0) = 1.09, 



while for p/p* = 4 the two expressions are within 1% of 
each other. If in addition we assume that T{p) « —£,{p) 
then the expression in Eq. (|3|) is again very similar to 
the form proposed by Joanny et a/Jl3 for two spheres (3 
being replaced by tt). Our arguments therefore provide 
an a-posteori justification of the validity of their poten- 
tial for the semi-dilute regime. In the dilute regime it will 
overestimate the attraction in the same way as Eq. jl^ ) 
does. 
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FIG. 7. Scaled depletion potential V"'-''\x) / (jiRc) for 
interacting polymers. From top to bottom the densities are 
p/p* = 0.2, 0.5, 1, 2, 5 and 10 respectively. In the semi-dilute 
regime the well-depth increases like p^-"^"^^ while the range 
decreases as p""'^^". Inset: Scaled well depth at contact 
V^'^'^^ (Q) / {-kRc) for the interacting polymers (solid line), com- 
pared to the s cale d well depth Fao (0)/(7r7?c) for the AO po- 
tential of Eq. (^ with Rao = Rao gi^^n by Eq. (0). The 



scaled AO potentials are shown for Rc = IQRg, (dotted Hne), 
Rc = 2Rg (dashed line) and Rc — 0.5Rg (long-dashed line). 
In contrast to the Derjaguin approximation expressions, the 
AO potential does not satisfy prefect scaling with Rc- The 
well-depth scales linearly in p. ;c owff 

independent of p. 



The range is 2RJq, and is 



Given that we now have a reasonably accurate deple- 
tion potential between two spheres, namely Eq. (p^, it 
seems fruitful to examine how this expression varies with 
density. We do this in Fig. ^ for a number of densities 
in the dilute and the semi-dilute regimes. Since within 
the Derjaguin approximation the depletion potential ( p^ ) 
always scales with Rc, the curves in the plot can be used 
for all size ratios. Of course one should keep in mind 
that the Derjaguin approximation becomes progressively 
less accurate for decreasing Rc/Rg. With this caveat in 
mind, the well-depth at contact for polymers in a good 
solvent scales as: 



vl(Zv-\) 



„0.770 



(14) 



in the semi-dilute regime. This is exactly half the expo- 
nent with which the well-depth at contact scales for two 
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walls in the semi-dilute regime (see Eq. (^), which can 
be understood as follows: Within the Derjaguin approxi- 
mation the force scales as F(0) = nRcW{0) ^2^/(3"-!); 
the potential is then obtained through Eq. by inte- 
grating this force over a range D{p) ~ p-'^/(3''-i) . These 
two effects result in the scaling seen in Eq. (^4|). On the 
other hand, by nature of the Derjaguin approximation, 
the range is the same as that shown in Fig. H for the 
wall- wall case, i.e. it decreases as D^{p) ~ p~^^'^^. 

An expression similar to Eq. (|l^) has,been derived by 
Fleer, Scheutjens and Vincent (FSV)ll3. They replace 
Rao in the full AO potential (|Ai1) by 2f(el and the 
density p by n(/3). In its original application^, the FSV 
theory was used with a self consistent field theory to cal- 
culate T[p), and a simple Flory Huggins prescription for 
n(/9). This leads to V{p) ~ /9~2 scaling behavior which 
is more appropriate for polymers near the theta pointtJ. 
However, one could easily generalize the FSV approach 
and include the r(p) and n(p) appropriate for polymers 
in a good solvent. For the semi-dilute regime the FSV ap- 
proach would then lead to an underestimate of the range 
by a factor 2/3, while for small g, where the Derjaguin 
approximation is quite accurate, it would underestimate 
the attraction at contact by a factor of about 4/9. On 
the other hand, in the dilute regime, this approach should 
perform better, although for increasing q it is expected to 
overestimate the attraction at contact in a very similar 
way to the effects seen for the AO model in Fig. |l^ (a). 
This is because the FSV model ignores the deformation 
of polymers near a spherical surface. It would be better 
to use the D{p) appropriate for a spherical particle; this 
results in a potential similar in spirit to an-jpteresting 
recent proposal by Tuinier and LekkerkerkeiO. 



D. An accurate semi-empirical depletion potential 

As can be seen in Fig. U, the Derjaguin approximation 
seems to be particularly accurate for the force at contact. 
In fact, when we plot F(0)/(7r i?c) v.s. p/ p* in Fig. ||, the 
results for the three size ratios are very close to each 
other at each of the densities studied, and their value 
is well described by 14^(0) — —2j^{p), as expected for 
the Derjaguin approximation (see Eq. O). Furthermore, 
we analyzedr-Uie, computer simulation data of Dickman 
and YethirajESH, which used the fluctuation bond model 
(FBM) for size ratios ranging from q = 1.3 to g = 3.5, 
and find similar agreement in Fig. ^. Keep in mind that 
the FBM results are for shorter chains (L = 20 — 100), 
and that the error bars appear to be larger than those in 
our simulations. Since for q = 3.5 the Derjaguin approx- 
imation (|l^) should not be valid at all, this does suggest 
that the expression 



F(0) = TrR,W{0) = -2^i?e7^(p), 



(15) 



may hold because of an (approximate) sum-rule, and that 
the agreement is not fortuitous. 
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FIG. 8. 

Comparison of the force at contact with the Der- 
jaguin expression from Eq. (11). The circles are taken 
from the SAW simulations depicted in Fig. ^. On the 
scale of this graph, the three size ratios lead to virtu- 
ally identical results at each density. The triangles de- 



note F(0)/(7r i?c) from the fluctuation bond model sim- 
ulations of Yethiraj and Dickmann3. The diamonds are 
from previous L = 100 SAW siraulations of W{0), the 
interaction between two platea23'E^. The curve denotes 
W{0) = — 27u,(p), the expected result for the Derjaguin 
approximation. 

The surprising accuracy of expression ( [l5| ) for the force 
at contact can be exploited to derive a more accurate 
effective potential than that found by the direct appli- 
cation of the Derjaguin approximation. We assume the 
same linear force between two spheres as found in the 
previous two sections, but now with a modified range 



Ds{p) = D^{p)- 



Ro 



(16) 



where the effective shrinking of the range due to the de- 
formation of the poly mer s around the HS has been taken 
into account via Eq. (A3). Although this deformation ef- 
fect is calculated for ideal polymers, we expect it to be a 
good first approximation for the relative deformation of 
interacting polymers around a sphere. The resulting sim- 
ple linear expression for the force can now be integrated 
up via Eq. (nffl) to give the depletion potential: 



(17) 



for X < Da{p). Vs{x) — for x > Ds(p). Our new semi- 
empirical potential works remarkably well, as can be seen 
in Fig. H for the three size ratios. Even for q = 1.68, 
corresponding to the smallest colloids, this potential is 
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significantly better than the Derjaguin expression (12), 
suggesting that Eq. (|l^) can be considered a nearly quan- 
titative representation of the depletion potential induced 
between two spheres by SAW polymers for a surprisingly 
wide range of q' s. Since the semi-empirical potential ( p7| ) 
reduces to the regular Derjaguin form (|l2[), if Ds{p) is 
replaced by D^{p), its scaling properties with density 
should be the same as those discussed in the previous 
subsections. 



E. Comparison with theories for non-interacting 
polymers 

It is interesting to compare the results for interacting 
polymers with the results for non- interacting polymers. 
Firstly, just as was found for two walls, the range for 
the depletion potential is independent of density for non- 
interacting-polymers. This is true both for the simpler 
AO models, as well as fntpjoce sophisticated theories 
and computer simulationsEJ'oS. For a given i?g, using 
a theory based on non-interacting polymers will always 
overestimate the range compared to the interacting case. 
Since for our accurate semi-empirical potential the range 
Ds {p) is related to D^, (p) by a density independent factor 
related to g, the relative overestimate in the range for the 
spherical case should be close to that found when com- 
paring a non-interacting theory to an interacting one for 
the range of the depletion potential between two plates, 
as done in Fig. ^. 

For the case of two walls the absolute value of the well- 
depth at contact given by ideal polymers or the AO model 
was underestimated w.r.t. the interacting case. Within 
the Derjaguin picture of Eq. (^l|), this implies that non- 
interacting polymer models will underestimate the force 
at contact -F(O) for two spheres in a bath of interacting 
polymers by a similar factor to that shown for W{Q) in 
Fig. 0. 

In the inset of Fig. ^ we compare the scaled well-depth 
at contact V^'^'^^{Q)/['kRc) for the Derjaguin approxima- 
tion with that of the AO model, given in Eq. (Al), where 
an effective range parameter R^/q^ given by Eq. (A2) in 



Appendix A, was used to take into account the defor- 
mation of the polymers around the colloids. For a given 
i?c, the well-depth Vao{^) scales linearly with p, which is 
now an overestimate compared to the result for interact- 
ing polymers. The origin of this behavior can be easily 
understood: while the strength of the force is strongly 
underestimated by the AO model, the range is overesti- 
mated. Therefore the integral in Eq. ( ]To| ) shows a com- 
pensation of errors. This explains why, at low densities, 
the AO approximation well depth is very close to the in- 
teracting polymer result. For example for q — 0.1 we 
find for p/p* = 1 that Vs{0) /Vao{'^) « 0.93 while for 
p/p* = 10 this ratio is still only 0.56. Attempting a 
naive improvement by replacing 11 (p) = p in the AO ex- 
pression by the true pressure of a polymer solution results 



in a worse approximation for all densities. 



F. Polymers as soft colloids 

When modeling polymer-colloid mixtures, the colloids 
are usually treated as single particles. It thus seems nat- 
ural to attempt the same for the polymers. We have 
recently ■succeeded in modeling linear polymers as "soft 
colloids"E3~cdi£3, where each polymer is replaced by a sin- 
gle particle centered on its CM. Different coils interact 
via an effective pair potential acting between their CM. 
These effective interactions v{r; p) between the CM of 
the polymers were extracted from the radial distribution 
functions g{r) with an Ornstein Zernike inversion tech- 
nique. This procedure is justified by a theorem which 
states that for any given g(r) and p there exists a single 
unique pair-potentj.al .which reproduces that radial dis- 
tribution functioiiE3'E3. Our w(r; p) have a near Gaussian 
shape, with an amplitude of order 2kBT and a range of 
the order of the radius of gyration Rg. For p/p* < 2 
accurate analytic forms as a function of r/Rg and p/p* 
are availablecd. 

Our input g{rys were generated by computer simula- 
tions of L = 500 SAW polymers on a cubic lattice, and 
so the density dependent effective potentials we derived 
will by definition reproduce the same CM pair structure 
as found in the underlying SAW polymer system. Wh^ 
used in conjunction with the compressibility equatioi 
these potentials also prpvidje a very good representation 
of the equation of stateEScil. 

A similar coarse-graining procedure was followed to de- 
scribe polymers near a flat wallE3Ej, where even ideal 
polymers form a depletion layer. If these were to be 
modeled as single particles, there would be no polymer- 
polymer interaction, but there would still be an effective 
polymer- wall interaction (j){z). We used direct simula- 
tions of L = 500 SAW polymers to extract the density 
profile p(z) near a hard wall. The 4>{zys which repro- 
duce the density profile at each given bulk density Pi W - er e 
calculated using a wall-Ornstein-Zernike procedureHE3. 
Again accurate analytic fits for densities p/p* < 2 are 
availablecil. 

Since our effective polymer-polymer potentials provide 
a very accurate representation of the pressure n(p), while 
the polymer-wall interactions are constrained to repro- 
duce the correct density profile p{z), and therefore the 
correct adsorption f (p), Eq. (||) implies that our soft- 
colloid approach will reproduce the correct wall-fluid sur- 
face tension 7^(p). This explains why our sofL-piUoid 
approach reproduces the correct value for iy(0)Ea'L3, the 
depletion potential at contact between two walls, since 
Eq. implies that this can be expressed completely 
in terms of 7u,(p). Similarly the slope of the potential 
is given by the osmotic pressure 11 (p) which is also ac- 
curately represented in the soft-colloid approach. De- 
viations w.r.t. direct simulations were only observed for 
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larger distances z, 
duced a more pro* 
tion potentials .EaE 



where the soft-colloid approach pro- 
punced repulsive bump in the deple- 



*^\\\ ' 

^\\\ 


'1 




q=1.68 







comes more pronounced with decreasing sphere size (this 
doesn't happen in the monomer representation of course) . 
Similarly, the effective potentials 4>{r), derived with our 
inverse OZ approach, become relatively softer with de- 
creasing Rc, as shown in Fig. Bpl. For a fixed q the 
potentials become more repulsive for increasing density, 
just as was found for the polymer-wall casa2j. Note that 
this behavior is the opposite to that of the density pro- 
files. Again, because these potentials are constrained to 
give the correct density profiles and related adsorptions, 
this soft colloid picture should also correctly reproduce 
the surface tension 7., (p) for a single sphere immersed in 
a bath of interacting polymers. The same should hold 
for the related one-body insertion free energies 
described in paper I. 




FIG. 9. The effective sphere-polymer potential (t>{r) as 
a function of r, the distance to the centre of the sphere for 
increasing polymer concentration (from left to right). These 
(j){r) are derived via an Ornstein-Zernike inversion technique 
such tljiat they reproduce the CM density profiles around each 
spherecil. The arrows denote the hard core radius of the 
spheres. 



To calculate the depiction potential between spheres, 
we need to first derive the effective sphere-polymer CM 
potential (f>{r) which would exactly reproduce the density 
profile p{r) around a single sphere of radius Rc- This 
can again be done using an Ornstein-Zernike (OZ) tech- 
nique to invert the density profilescj. These density pro- 
files were calculated in paper I, where we found that the 
range of the profiles shrinks with increasing density. The 
CM profiles can penetrate into the spheres because the 
polymers can deform around them, an effect which be- 



1. Direct simulations of depletion potentials for polymers as 
soft colloids 

The simulation of the depletion potentials induced 
by the "soft colloids" proceeds in a similar manner to 
the simu lations with SAW polymers, described in see- 
However, the implementation is much simpler. 



tion III A 



because the interactions are smooth and both types of 
particles are in continuous space. The simulations are 
also about two orders of magnitude faster, since each 
polymer is represented by only one particle, instead of 
the 500 units per polymer used for the SAW model. 
These simulations are compared for q = 1.05 in Fig. |l^ 
to the direct SAW simulations. The agreement is very 
good for p/p* — 0.58 and reasonable for p/p* = 1.16. 
However for the highest density, p/p* — 2.32, the soft- 
colloid approach breaks down. Exactly the same trend 
was found for the other two size ratios (not shown here). 
Since the soft-colloid approach reproduces the correct 
one-body densities p{r), and the correct one-body free 
energy F"'* (p) of immersing a single sphere into a bath 
of polymers for all densities p/p* in the dilute and semi- 
dilute regimes, it is perhaps surprising that for the case 
of two spheres, it breaks down rather abruptly for higher 
densities. To further investigate this issue, we plot the 
polymer CM density profiles for two spheres fixed at dis- 
tance r = 2.6Rg for q = 1.05 and p/p* = 2.32 in Fig. |l|. 
This corresponds to a distance where the SAW V{r) has 
effectively gone to zero. At this short distance the SAW 
polymers clearly penetrate more easily in between the 
two colloids than the soft particles do, resulting to less at- 
traction between the colloids. A,similar effect was found 
for polymers between two wallsEj, but there the fact that 
sum-rules constrain the value of the potential at contact 
means that this difference manifested itself instead in an 
enhanced repulsive bump in the potential. In both cases, 
we attribute the error in the soft-colloid approach at 
strong confinement and higher densities to a hceakdown 
in the "potential superposition approximation"!^, i.e. the 
assumption that a polymer confined between two surfaces 
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feels an interaction which is simply the sum of the two 
interactions (/)(r) with each separate surface. While this 
would be correct for simple fluids, it is not correct here, 
since the close presence of one wall affects the interaction 
of the polymers with the second wall. This is essentially 
due to the deformability of the polymers, which is not 
correctly taken into account for strong confinement. 




SAW simulations p=0.58 
SAW simulations p=1.16 
SAW simulations p=2.32 
soft coiioids p=0.58 
soft coiioids p=1.16 
soft coiioids p=2.32 
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x/R„ 



1.5 



FIG. 10. Comparison of the direct SAW simulations of the 
depletion potential for q — 1.05 to direct simulations within 
the polymers as soft colloids approach. 




8cm 



FIG. 11. Comparison of the direct L = 500 SAW 
simulations of the polymer CM density profile to direct 
simulations within the polymers as soft colloids approach 
for q = 1.05 and p/ p* = 2.32. One sphere is at 0, and the 
other is at r = 2.6Rg. The four graphs show the normal- 
ized density at a distance r, averaged over (0 — lOdeg), 
(10 - 20deg), (20 - 30deg), and (30 - 40deg) w.r.t. the 
particle at the origin. The SAW polymers penetrate more 
easily between the colloids than the effective soft-colloids 
do. 



When we do a similar analysis of the densities as done 
in Fig. ^ for p/p* = 1.16 however, the differences in 
the density profiles are very small, suggesting that the 
soft-particle picture is not yet breaking down. Coarse- 
graining polymers as soft colloids is expected to be most 
useful for q < 1, since once the polymers become much 
larger than the colloids, they can easily wrap around the 
colloids, an effect not easily treated in our coarse-graining 
scheme. Since for q < I, we have found that phase- 
separation sets in below p/p* « lE3, the limit of physical 
stability is reached before we encounter problems with 
our coarse-graining scheme for polymer-colloid mixtures. 



2. Depletion potentials from the superposition approximation 

To understand further the nature of the depletion po- 
tentials between the colloidal particles, mediated by the 
polymers in the "soft colloid" picture, we now turn to the 
superposition approximation, which expresses the two 
body depletion interactions in terms of one-body proper- 
ties of a single colloid in a bath of soft polymers. 

If one colloidal particle is fixed at the origin, and one 
is fixed at r, then one can define an inhomogeneous one- 
particle density of the soft colloids p'""^' (r'; 0, r) for that 
fixed configuration of the two colloidal particles. The ef- 
fective force induced between twapGolloidal HS by the soft 
colloids can then be written ascBE^I: 



(18) 



were (f)(r) is the interaction between the HS particles and 
the soft colloids. This expression has an obvious intu- 
itive interpretation. The total force on a HS particle at 
is the average of the sum of its interactions with all the 
soft colloids. If there are no other HS particles, then the 
equilibrium distribution will be spherically symmetric, 
and this force will be zero. However, if there is another 
HS particle at r, it will perturb the density distribution 
of the soft colloids, which will in turn result in a different 
total force acting on the particle at 0. By definition, this 
is the effective or depletion force F(r) acting on a particle 
at 0, induced by the presence of a particle at r. 

Eq. ( p^ ) is in principle exact, but requires knowledge of 
the soft colloid density around two colloidal HS. By ap- 
proximating the full (cylindrically symmetric) one-body 
density p(^^(r';0,r) by a superposition of the,one-body 
density p{r) around an isolated single sphereE^: 



p«(r';0,r) = p(r')p(|r-r'|)/p. 



(19) 



the effective depletion force (|8|) can be entirely expressed 
in terms of the (radially symmetric) problem of a single 
colloidal sphere immersed in a polymer solution. The 
results of such a calculation for V(r) are compared in 
Fig. |l^ for q = 1.05. For all three densities, the super- 
position approximation closely follows the direct simu- 
lations for the soft-colloids. This implies that the full 
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one-body density /5^^^(r'; 0, r) does not differ much from 
the simple superposition of Eq. (|l9|), even for the high- 
est density considered, p/ = 2.32. The reason for this 
is most hkely that the effective polymer-polymer CM 
interaction v(r) is rather weak. In a recent study, a 
similar good performance for the superposition approx- 
imation was found for a low density hard-core deple- 
tant wlien the HS-depletant interaction was fairly long 
rangedcj. Good agreement was also found in a similar 
study of star-polymer colloid|anixturesEil where accurate 
star-polymer-wall potent4alsE3 are derived based on a 
coarse- graining approachS In contrast, if the density 
of a hard-core depletant increases, then the superposi- 
tion approximatioja is known to break down for increasing 
packing fractionE^O. 



teractions are therefore quite sensitive to the exact form 
of 0(r). 

In conclusion then, the superposition approximation 
works remarkably well for p/p* < 1, the regime where 
the soft-colloid picture provides the most accurate repre- 
sentation of the true depletion potentials. Since the den- 
sity profiles p{r) can be easily calculated by our wall-OZ 
approach, this means that one only needs the polymer- 
colloid interaction (j)(r) and the polymer-polymer interac- 
tion v{r) as input to reliably calculate the full two-body 
depletion interactions from the superposition approxima- 
tion. 



IV. SECOND VIRIAL COEFFICIENTS 
PHASE-DIAGRAMS 








soft colloids p=0.58 

- soft colloids p=1.16 

- soft colloids p=2.32 
superposition p=0.58 
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FIG. 12. Comparison of simulations of the depletion po- 
tential for q = 1.05 within the soft-colloid approach to the 
superposition approximation of Eqs. (^) and ^). The thin 
solid line denotes the superposition calculation with p{r ) from 
p/p* = 2.32, but the potential (/)(r) from the p/p* — 1.16 
calculation. This illustrates the sensitivity of the depletion 
potentials to colloid-polymer CM interaction (j>{r). 



The fact that the superposition approximation works 
well for the soft colloids at the higher densities, also raises 
an interesting question as to why the soft-colloid picture 
itself breaks down at the higher densities. By construc- 
tion, the one-body density profiles are identical for the 
full polymer case and the soft-colloid picture. Presum- 
ably the reason for the difference in the two-body deple- 
tion interactions can be traced either to a large deviation 
from superposition for the full polymer case, or else to a 
breakdown of the potential superposition approximation, 
as discussed in ref.tj. The influence of changing (j>{r) is 
illustrated in Fig. ^ where we used the (/)(r) appropriate 
for p/p* = 1.16 to calculate the depletion potential with 
the density profiles for p/p* — 2.32. This 4>{r) is less 
repulsive than the true (/'(r), and leads to an important 
difference in the depletion potential. These depletion in- 



A useful measure of the strength of an interaction is 
the second osmotic virial coefficient B2, since this is ex- 
perimentally accessible through ^measurement of the 
osmotic pressure at low densitiesEj'Ej. Such virial coeffi- 
cients are very sensitive to the strengthr^d nature of the 
depletion interactioncS. A recent studyO has shown that 
for many simple systems consisting of a HS like repul- 
sion with an additional attractive potential, the location 
of the liquid-gas or fluid-fluid critical point can be corre- 
lated with the point where the reduced virial coefficient 
= B2/B§^ w -1.5, where B^^ = IQttRI/?, is the 
virial coefficient of a HS system with radius Re- For our 
depletion systems, the reduced virial coefficient is given 
by 

3 



B*2 



1 



4i?3 



(exp[-y(a;)] - 1) , 



(20) 



The direct simulation results depicted in Fig ^ can be 
used to calculate B2 , and the results are shown in Fig |l^ 
for p/p* = 0.58, p/p* = 1.16, and for p/p* = 2.32. As 
expected, for a given density p/p*, the B2 become pro- 
gressively more negative with decreasing q. Although we 
don't include any explicit error bars in our plots, these 
may be rather large for increasing attraction because they 
appear exponentially in Eq. (pO|). 

We found in the previous sections that a simple semi- 
empirical potential Vs{x), given by Eq. (17), provided a 
nearly quantitative description of the depletion potentials 
between two spheres. It turns out that for this potential 
the virial coefficient can be integrated analytically: 



B2 

T^HS 
"2 



1 



X^TiR^J'^W'^l 



(2^/RcW {~6D{D + 2Rc 



'SD cxp 



-RcWD 
2 



{D + AR,) + t:R^W{D + 2R^)^R, 



t:RcW{D + 2Rcf)Erfi 



TTRrDW 



where W = W{x = 0) and D = Ds{p), i.e. the ex- 
plicit density dependence has been suppressed for nota- 
tional clarity. The imaginary error function Erfi[z] is 
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defined as Erf[iz]/i; its value is real. As can be seen 
in Fig. |l^, this expression compares well with the direct 
simulation results, demonstrating that the potential ( |l7| ) 
is a good approximation to the true depletion potential. 
For this reason, we also plot the theoretical B2 for two 
smaller q's. The density p/p* at which one would expect 
a (metastable) critical point and a related fluid-fluid de- 
mixing decreases with decreasing q. To compare with 
experiments, one should keep in mind that our parame- 
ter p/p* would be equal to the density of a reservoir of 
polymers kept at the saHie chemical potential as a full 
polymer-colloid mixturai3. 



GO 
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FIG. 13. Virial coefficients; the symbols represent results 
from the computer simulations of L=500 SAW polymers for 
q — 0.67 (circles), q = 1.05 (squares), and q = 1.68 (trian- 
gles). The straight lines are from Eq. (pl[). The good agree- 
ment shows that Eq. (|l^) provides a good representation of 
the depletion potential. 

If the pair-potentials are of sufficiently short range, the 
fluid-fluid transition becomes metastable-sv.r.t. the fluid- 
solid transition of the larger particlesQEiH, which widens 
to a much larger. It has been recently found that this 
generic widening of the fluid-solid transition is priipaji- 
ily driven by the value of the potential at contactEScI; 
other details such as the shape and range of the attrac- 
tive part of the potential are less important. Roughly 
speaking, the fluid-solid liquidus curve widens to about 
half the packing fraction of the hard-sphere freezing tran- 
sition when V{0) « 2.5 ± 0.50. When combined with the 
criterion that B2 ~ —1.5 at the critical density, this can 
be used for a rough prediction of the q below which the 
fluid-fluid transition is metastable w.r.t. the fluid-solid 
one. The minimum q for which the fluid-fluid transition 
is still stable is betwee]a.|0p3| |a«d 0.5, which is ccwpstcnt 
with other calculationsLfoEJ'tll, and experimentsEillj. Of 
course one should keep in mind that the larger the size 
ratio q, the more important three-body and higher or- 
der interactions becomeE^IIj. These pair-potential crite- 
ria should therefore become less accurate for increasing 

q- 



Both P4f~l IheoryES and integral equation 
approachesOO predict that for large enough q, B2 
should go through a minimum as a function of p/p*, 
although the two theories differ signiflcantly on quantita- 
tive details. Experiments seem to show similar effectsES. 
Remarkably, expression (|2^), which is based on the sim- 
ple potential of Eq. (^^ also shows a minimum, as shown 
in Fig. |lj. Even more surprisingly, considering how 
the potential was derived, is that the minimum is near 
p/p* — 1, as predicted by RG theory, and that it closely 
follows the asymptotiC|-4|aw min{B2) = 1 — 0.5g°'^°^, 
derived by EisenrieglerEil (see the inset of Fig. |l4). In 
contrast, when the simpler Derjaguin expression is 
used for _B|, no minimum as a function of p/p* appears 
for large q. At present it is not clear why our expression 
for B2 , Eq. ( ^ ) , should agree so well with the RG results 
in this large Rg/Rc regime. 
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FIG. 14. Virial coefficient from Eq. ( [ill ) ^ ^ function of 
polymer density for large size ratios; q = 5 (solid line), q = IQ 
(dotted line) q = 20 (dashed line), g = 50 (long-dashed line) 
and q — 100 (dot-dashed line). Inset; Comparison of the min- 
imum of B2 as a function of q to the asymptotic RG result 
min{B2) = 1 — 0.5g'' ''''^ . The agreement is remarkable. 



V. DISCUSSION AND CONCLUSION 



We have used a combination of computer simulations 
and theory to derive the depletion interactions induced 
by excluded volume polymers. Our main computer sim- 
ulation results are depicted in Fig. |4| and Fig. ^ and our 
main theoretical results are the depletion potential for 
two plates, given by Eq. (|^) and the depletion potential 
for two spheres, given by Eq. ([TtI). 

For two plates, we found that a simple linear depletion 
potential (|) with a well-depth of W{0) = 2'yw{p), and 
a slope of n(/9) is very accurate. The range decreases 
with density, with the largest relative change occurring 
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in the dilute regime. In the semi-dilute regime, the well- Eq. ( pl| ) works so well in this limit. 



depth scales as W{0) ~ p , while the range scales as 



-0.770 



The depletion potential also takes on a 



simpler functional form, given by Eq. (|6|). These results 
differ significantly from theories based on ideal polymers. 

For two spheres, we find that the simple Derjaguin ap- 
proximation works much better than expected due to a 
cancellation of errors related to the deformation of the 
polymers around a sphere. Interestingly, the Derjaguin 
expression for the force at contact seems particularly ac- 
curate, even for very large size ratios q. We used this 
observation, together with a correction for the Derjaguin 
range, to derive the depletion potential Vs{r) (Eq. (p7[)), 
which was shown to agree quantitatively with our di- 
rect computer simulation results. While the excellent 
agreement with simulations is gratifying, further study 
is needed to understand why our main ansatz, that the 
force at contact is accurately approximated by the Der- 
jaguin approximation, works so well for such a wide range 
of q. 

The range of our potential scales with density in the 
same way that the range for the two plates does. But in 
contrast to the case of two plates, where the well depth at 
contact is more attractive than that for ideal polymers, 
the well depth at contact for two spheres is less attractive 
than that predicted by ideal polymer theories. For ex- 
ample, in the semi-dilute regime, Vs{Q) ~ —p^ ™, com- 
pared to the scaling of the AO potential: Vao(O) ~ —p. 
In other words, ideal polymer theories overestimate both 
the range and the well-depth of the depletion potential 
between two spheres, induced by interacting polymers. 

We also calculated the depletion potential between 
spheres induced by polymers coarse-grained as "soft col- 
loids". For the dilute regime, good agreement with the 
direct simulations of SAW polymers was achieved. How- 
ever, for the highest density, p/p* — 2.32, the polymers as 
soft colloids approach strongly overestimates the depth of 
the depletion potential. We traced this discrepancy to a 
breakdown of the potential superposition approximation, 
something also found for two platesE3. However, since 
our coarse-graining scheme is mainly useful for q < 1, 
where phase-separation sets in for p/p* < 1, this break- 
down at high density is not relevant to the equilibrium 
phase behavior of colloid-polymer systems. In fact, we 
have recently used this scheme to calculate the phase- 
separation binodals ©f, polymer-colloid mixtures for sev- 
eral size ratio's q < llI3. Direct simulations of SAW poly- 
mers and colloids would be significantly more expensive. 
The advantage over directly using a pair-potential such 
as Eq. is that many- body interactions are also ef- 

fectively taken into accountlla. 

Virial coefficients can provide a sensitive measure of 
the quality of an effective potential. We were able to in- 
tegrate Eq. (|l^) to derive an analytical representation of 
the virial coefficient which compares very well with 
the direct simulation results. Surprisingly, this same an- 
alytic form agrees quantitatively with RG results in the 
limit g 3> 1. Further study is needed to clarify why 



A careful comparison with experiments is complicated 
by the fact that most measurements are of phase be- 
havior or structure, which depend only indirectly on the 
depletion interactions. The most extensive direct mea- 
surements of the depletion .pptential between spheres, 
carried out by Verma et a/Jl3, were on a system with 
semi-flexible polymers (DNA) which do not follow the 
same scaling laws as the SAW polymers we treated in 
this paper. However, given a good model for the equa- 
tion of state and the adsorption r(p) for these polymers, 
it would be straightforward to calculate the depletion po- 
tentials using methods derived in this paper. A detailed 
comparison with-jepcperiments will be the subject of an 
upcoming reviewO. 

This paper has concerned itself with the case of in- 
teracting, non-adsorbing polymers in a good solvent. By 
taking Eqs. (|l|), (|), (||), and ^ together, if follows that 
accurate depletion potentials can be calculated with only 
a knowledge of the adsorption T{p) and the equation of 
state H(p). This implies that the effect of changing sol- 
vent quality, or of adding attractions or repulsions be- 
tween the colloids and the polymers can be understood 
by seeing how they effect these two quantities. Note 
that this is different from binary colloid mixtures, where 
the effects of added interparticle attractions and repul- 
sions have_a more subtle effect on the resultant depletion 
potentialsCJ. 

It is not hard to see that for polymers in a good solvent, 
adding an attractive (repulsive) polymer-colloid interac- 
tion results in a less (more) attractive depletion potential, 
and possibly even to repulsive effective interactions. This 
follows since H(p) is unchanged, and only T{p) is affected 
in an obvious way. An example where this may be rele- 
vant concerns the common practice to stabilize colloidal 
particles by a short polymer brush. The interaction of 
the free, polymers with the polymer brush can be quite 
subtleL3, but our theory implies that it is only the over- 
all change in the adsorption which effects the depletion 
potential, and by extension, the phase behavior. 

Changing the solvent quality has a less obvious effect, 
since both H(p) and Tp) are affected. For a given density, 
we expect the absolute magnitudes of both to decrease, 
so that the strength of the depletion potential should de- 
crease as well. Of course, for strong deviations from the 
case of non-adsorbing polymers in a good solvent, some 
of the simplifying assumptions that went in to the deriva- 
tion of the depletion potentials may begin to break down. 
For example, for very strong, adsorption, the system may 
show bridging flocculationcB, or gelationO. A rich phe- 
nomenology may be expected; this will be the subject of 
an upcoming investigation^. 
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APPENDIX A: DERJAGUIN APPROXIMATION 
FOR THE ASAKURA-OOSAWA MODEL 



two walls, one should take Rao = {'^/ V^)Rg, which cor- 
responds to approximating the depletion layer of ideal 
polymers by a step-function with the same depletion vol- 
ume. This then results in a well-depth which is equal to 
that found for ideal polymeraSa. 

However, for polymers of a given size Rg, the width of 
the depletion layer around a single sphere will decrease 
with decreasing Rc, because the pelft'saers can partially 
wrap around the colloidal particleg23C3. By equating the 
one-body insertion free energy for a HS sphere into a 
bath of AO model particles to that found for inserting a 
HS into a bath of ideal polymers, we found, in paper I, 
an exact analytic prescription for the effective parameter 



1. Comparison for fixed parameter Rao 



The Asakura Oosawa model, where the ideal polymers 
are modeled as interpenetra|ble spheres of radius Rao, 
was first introduced in 1958EI. Between two plates this 
results in the linear depletion potential similar to Eq. (|^), 
with a range Rao, while between two spheres it results 
in the well-known Asakura-Oosawa form 



R 



eff 
AO 



Rn 




(A3) 



VAo{x) = -p-{ac,) 



+ 



1 /x + 2Rc 
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< X < 2Rao (A1) 



which is exact within the confines of the AO model. Here 
we defined acp = Rc + Rao- Similarly, the use of the Der- 
jaguin approximation ( pi] ) together with the linear AO 
potential between two plates results in: 



V^S'^'{x) = -p^Rc{2R 



Ao-xy 0<x<2Rao- 

(A2) 



Comparison of the two 
Vao{0) = -p2tt{R 



r2 
^AO 



expressions 
2i?io/3) 



at contact, where 
and V^^''\Ol 



—p2ttRcR\q, shows that the Derjaguin expression ( [A2| ) 
underestimates the well depth at contact w.r.t. the exact 
expression( Al), a discrepancy which gets relatively worse 
for decreasing R^/ Rao ■ This is illustrated in Fig [T5|(a) 
for the same size ratios we used for the interacting poly- 
mers. Since the Derjaguin approximation is only valid 
for large Rc/Rao, this breakdown is expected. 



2. Comparison for effective parameter R^jIq {Rg , Rc) 



Before using the AO model to describe polymers, one 
needs to fix the parameter Rao- rA common prescrip- 
tion has been to take Rao ~ RM- However, between 



which decreases with decreasing Rc as expected. Just be- 
cause the one-body terms are equal does not mean that 
the same prescription for Rao will work fnrJthe two-body 
terms. But direct computer simulationscj have shown 



that using R^'/l in the full AO expression (|A1[) results in 



a much better approximation of the depletion potential 
between two spheres in a bath of ideal polymers, than 
using a prescription which fixes Rao independently of 
Rc- 

With this in mind, we now compare again the Der- 
jaguin expressi on (|A2| ) with Rao = S/y^/Jg to the full 
AO expression (|AlD but now with an effective AO radius 
R'^Iq from Eq. (|A3[) appropriate to the desired Rc. As 
demonstrated in Fig. |l^(b), the Derjaguin approach is 
now a more faithful approximation of the full AO model 
depletion potential, which is itself a better approximation 
of the true depletion potential induced by ideal polymers. 
Although the Derjaguin approach doesn't capture the 
slight reduction in range, the well-depth is much closer 
to the AO value. 



Of course this improvement arises from a cancellation 
of errors. However, for the interacting polymers the de- 
pletion potential between two walls is also well described 
by a linear form, and we expect the same deformation ef- 
fect to be relevant for polymers near spheres. Therefore, 
we might hope for a similar cancellation of errors when 
the Derjaguin expression is used for interacting polymers. 
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(a) 



V»o(x) 
C(x) 




x/R, 

FIG. 15. (a) Com parison of the full Asakura-Oosawa 
model interaction (Al) wi th R ao = Z/y'Tr to the Derjaguin 



approximation expression (A2) with the same Rao- The hard 
sphere sizes are the same as in Fig. ^. The potentials, from 
bottom to top, correspond to q — 1.68, q = 1.05 and q — 0.67 
respectively. The polymer density is p/p* ~ 1. (b) The 
same comparison as in (a), but now the decrease of the de- 
pletion layer range with decreasing sphere size is taken into 
account by using the effective AO parameter R'^q, defined in 



Eq. (^), in the full AO expression (Al). For the Derjaguin 
approximation expression Rao = is unchanged. The 

agreement is now much better. 
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